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STABILITY OF SAS-A DUAL SPIN SPACECRAFT 
WITH ENERGY DISSIPATION ON THE MOMENTUM WHEEL 

Peter M. Bainum 


ABSTRACT 

The attitude stability of the SAS-A satellite with damping in the momentum wheel 
as well as the "despun** portion is analyzed. Wheel energy dissipation is modeled 
by assuming the wheel can flex with two degrees of freedom relative to the hub. 
The nonlinear attitude equations are derived for small wheel flexural motion and 
are a ninth order nonautonomous set. If the main body damper mass and wheel 
transverse moment of inertia are assumed small when compared with main sat- 
ellite masses and inertias, an averaging process can be used to determine the 
zeroth and first order secular perturbations on the behavior of the system mira- 
tion angle. From this a general analytic stability' criterion is established. A 
numerical evaluation of this criterion using SAS-A parameters and measured 
wheel damping data indicates that stability about a zero degree nutation angle is 
insured by a factor of 128 under normal operating conditions. 
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NOMENCLATURE 


A, B, C = main body moments of inertia about the x, y, z axes respectively 

A, B, C = composite moments of inei'tia about the x,y , z axes respectively 

b^ unit vectors along the x, y , z axes respectively (i = 1, 2, 3) 

b.' » unit vectors fixed to the nominal plane of the undetected wheel and 
rotating with it 

Cj = coefficients occurring in the steady state soliuions for <t>^, a,, 
dj , dj = the inner and outer radius of the rotor rim, respectively 
H - magnitude of system angular momentum vector H 
h = thiclmess of rotor disc 
= moment of inertia of main body about the b. axis 

I_ = moment of inertia of rotor about its b j axis 

= moment of inertia of the pendulous damper about the b^ axis 

” i 

K = the restoring spring constant of the torsion wire support 

= the restoring spring constant of the rotor 

K = K/mr^^ dimensionless form of K 

k = the nutation damping (rate) constant 

k_ == rotor damping constant 
1 \ 

^ = height of damper plane above x, z plane 
L. = the applied external torques about the b^ axis 
M = the mass of the main satellite and the rotor 
M = the total system mass 
m = the pendulum end mass 


vii 


nijj = the mass of the rotor 


q “ B/A 



r = radial coordinate of a differential mass on the rotor 
= the distance from the nominal spin (y ) axis to the pendulum hinge point 
Tj - the length of the pendulum 
s = spin rate of rotor relative to main body 

^ = critical value of s v/hen only one stable equilibrium angle exists 
T = kinetic energy 
t = time 

V = potential energy 

X, y, z = principal ares of main satellite 

y^ = displacement coordinate of from plane of undcflected rotor disc 
a^, a^= rotor deflection angles about the b^',E^' axes respectively 


r = 

y = the nutation angle, i.e. the angle between the b^ axis and H 

\ - [(B- A)n*I«j s]/a, in the zeroth order solution, frequency with which 
transverse angular velocity component rotate about b^ axis 

Q = nominal main bod / spin rate 

= in the zeroth order solution, the value of at some reference time 

cu. = angular velocities about the x, y, z axes respectively (i = 1, 2, 3) 

0 = position angle of the projection of dm^^ on the plane of the undeflected 
rotor disc 
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- nutation dampiir displacomont angle 
.0 = density of rotor disc 

T =: time constant ascocia:ed with nutation angle decay 
Superscript 

( )' refers to component in b,' axe® system 
( )’ Indicates differentiation with respect to time 

Subscript 

ave refers to averaged quanti ty 
b. , i refers to particular main body axes (i = 1, 2, 3) 


0 refers to initial state or cnuilibrium value 


STABIUTY OF SAS-A DUAL SPIN SPACECRAFT 
WITH ENERGY DISSIPATION ON THE MOMENTUM WHEEI 


I. INTRODUCTION 

NASA Goddard Space Flight Center Is currently directing the design and develop- 
ment of a dual-spin Small Astronomy Satellite (SAS-A) to be launched in the Fall 
1970. The satellite will have the capability of scanning the entire celestial sphere 
to determine the relative position of X-ray emitting sources with respect to the 
fixed position of the stars. It is important that the attituae cf the satellite be 
precisely known and maintained in order to accurately determine the location of 
the X-ray emitting sources. 

The damping and attitude stability of the SAS-A satellite with damping only on 
the slowly spinning main part was reported previously.^ The resulting differential 
equations of rotational motion when linearized were an autonomous set of fifth 
order equations with constant coefficients. Analytical stability criteria were de- 
veloped from these equations using the method of Routh-Hurwitz. 

Subsequent to the previous analysis, it has been demonstrated by static and dy- 
namic tests of the SAS-A momentum wheel that there was some energy dissipation 
in the shaft- momentum wheel assembly. The purpose of the present investiga- 
tion is to incorporate the effects of momentum wheel damping into the rotational 
equations of motion for the SAS-A spacecraft, and to analytically Investigate the 
attitude stability of such a system. A possible way to model the wheel energy 
dissipiition is to add a spring-mass-dashpot damper to the wheel similar to the 
model of Mingori.2 It is then necessary to relate the dissipation of the dashpot 
damper to that of the shaft- momentum wheel assembly. Another approach is 
that used by Sen^ and Fleisher^ in which the wheel is assumed to flex with t^^*o 
degrees of freedom with respect to the hub. 'Fhis latter alternative was selected 
as a more appropriate model since the results of the momentum w'heel tests 
could be directly used in determining appropriate wheel restoring and damping 
coefficients 

The present analysis is based upon the dual spin attitude control system con- 
sisting of a slowly rotating main part to which is attached the X-ray sensors, a 
high speed momentum wheel (rotor) whose spin axis is nominally parallel to the 
main body spin axis, and a pendulous type nutation damper attached to the main 
part and constrained to move in a plane which is perpendicular to the main body 
nominal spin axis. 
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:i. ANALYSIS 


The satellite, wheel, and damping systems are illustrated in Figure 1. The ele- 
ments of the attitude control system are: 1) the main part of the spacecraft, es- 
sentially a right circular cylinder w'here the nominal spin axis is the bj body 
axis, 2) a small momentum wheel or rotor assumed to be connected to the primary 
part near its center of mass, and whose spin axis is nominally parallel to the b^ 
body axis, and 3) a i^endulous-type nutation damper which is attocl^d to the main 
part and constrained to move in a plane a distance ■I above tl^ plane (deter- 

mined by the body axes perpendicular to the symmetry axis, b^). The dam|)er is 
hinged about a torsion wire support which offers a restoring torque in addition 


y SPIN AXIS 



FigL"’e 1. Eiements of SAS-A Attitude Control System 
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to the dissipative torque associated with the damper motion. {A recent paper 
describes the design objectives and functional description of this damping system."*) 
Ill addition the rotor is assumed to have two degrees of flexural freedom with 
respe :t to the hub. In other words the plane cf the rotor disc is allowed to deviate 
through small angular deflections from its nominal orientation which is pc:T*allel 
to the X , z plane. 

A. Equations of Motion 


The development of the equations of motion follows that of Reference 1, except 
now it is necessary to reformulate the expression for the kinetic energy of the 
rotor. The rotor is assumed to be spinning with a constant relative angular ve- 
locity magnitude, s, with respect to the main spacecraft. The nominal direction of 
♦he rotor spin axis is the direction of the b^ unit vector. 

For the development of the rotor kinetic energy consider an incremental mass, 
dmp , on the rotor disc which has a displacement coordinate from the nominal 
rotor plane (Figure 2). The b^', b^', b^' unit vectors are fixed to the nominal 
plane of the wheel and rotating with it. It has been assumed that the rotor is 
attached at or very near the center of mass of the main spacecraft -damper sys- 
tem. Neglecting any offset between the attachment point and the system mass 
center, the |>osltion vector of dm„ relative to the c.m. can be expressed for very 
small y displacement by: 


r - r sin0bj' + y b^ ^ r cos 4 b^' 


( 1 ) 


where 4 is the position angle of the projection of diOj^ on the plane of the unde- 
flected wheel disc and does not vary with time. Utilizing the relationship between 
the time rates of change of a vector in an inertial and rotating system, the inertial 
velocity of dm^^ relative to 0 can bo expressed 

Vq = r cos<^) b^' 

sin0) b^ (2) 

where are the inertial angular velocity components of the rotor. 

Sen and Fleisherd proceed to diiferontiate Equation (2) again with respect to time 
and obtain, from force and moment considerations, the two equations of motion 
for a wheel having two degr? of flexural freedom with respect to the hub under 
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WHEEL SPIN AXIS 


2 / t2 



PLANE OF UNDEFLECTED WHEEL 


Figure 2. Momenttm Wheel Deflection Geometry 


suitable assumptions of small deflections. Because it is of interest here to de- 
termine the complete rotational equations of motion for the system, wo will pro 
ceed to develop an expression for the kinetic energy of the rotor. 

The kinetic energy of the incremental mass dm may be expressed as: 


K 



If it is assumed that the rotor has two degrees of flexural freedom with respect 
to the hub, i.e. deflection about the b^' axis and deflection about b^', then 
the displacement coordinate can be written in terms of the flexure angles as: 


y^ = r sin0a^ •• rcos0a^ (4) 

for small deflections when the arc length can be approximated by the tangent. 

(For the general case of large deflections the position of dnij^ would have to be 
expressed in terms of spherical coordinates or an Euler angle sequence.) The 
differential mass can be related to the unit disc thickness, h, and uniform density, 
p, according to: 

dnij^ = p hr dr 6(p (5) 

After substituting Equations (2), (4), ana (5) into Equation (3) and integrating 
over the wheel volume. 



and m„ is the mass of the rotor rim; d., d are the inner and outer radius of the 

K * A 

rim, respectively. 



By writing the oo! in terms of the (see Figure 2), Equation (6) may be ex- 
panded in terms of the spacecraft body angular rates to yield; 

Trotor " 1 [iRj K* Ir, 

2 cos St +cc)j sin st) ~ 2 (- co^ sin st + co^ cos st] 

” 2 (“ sin st + o)^ cos stj sj 

“ 2 cos st + aij sin st) + s) 

^ “3') ^ “x^l 



In the limiting case, if the rotor is assumed to be completely rigid without 
flexural freedom with respect to the hub, a^=a^=a^ = a^ = 0, and. 


( 8 ) 


ROTOR 


2 


Jrj K* ^ Ir, W*"3') 


Equation (8) is seen to be identical to Equation (1) of Reference 1. 

An additional centrifugal force potential exists which is associated with con- 
servative forces tending to restore the deflected rotor rim to the nominal b^, b^ 
plane. With the aid oi Figure 2 it can be seen that the centrifugal force acting 
on dnij^ due to the rotation about the bj axis can be approximated, for a « 1 and 
y « r , by 


Fc = dm^ r(co2+ s)2 


( 9 ) 
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The component of tending to accelerate toward the plane of the undeflected 
wheel is: 


F 


c 


a 


F sina ^ 

C 






( 10 ) 


But for small deflections ra ^ so that 


Fc ^ r(a^sinc^+a^ cos0) dm^^ (11) 

To obtain the potential, consider the work done by F^ through the differential 
distance dy^^ cos a ^ dy^ = r^sin0da^+ cos j, “ 


dW = F^ • dy_ (11) 

After performing the multiple integration, and recalling the relationship that 
Irj ” "’r ^ ^ 2 )/^* following expression is obtained: 


ROTOR. 



* s)2 (a 2 


( 12 ) 


If it is now assumed that the wheel dissipative forces vary linearly with the 
flexural angular rates and<i^ , these forces can be derived from a Rayleigh 
dissipation function, 3, similar to that given in Reference 1, which has the form: 


3 



(13) 


where k^^is the equivalent linear viscous damping constant of the rotor. Wheel 
damping tests performed at the Applied Physics Laboratory verify that, to first 
order, there is a linear relationship between wheel dissipative forces and the 
flexural angular rates. (Likins et al.^ have recently discussed some of the dif- 
ficulties which can result due to damping nonlinearities in dual spin systems and 
have also presented a technique of relating an equivalent linear viscous damping 
coefficient to the actual nonlinear damping.) 
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A’ 


In addition to the centrifugal restoring forces on the rotor, there is also a re- 
storing effect due to the stiffness of the rotor material itself. For small de- 
flections and homogeneous wheel material these forces can be obtained from the 
following potential energy expression: 


p 





'•:*s 


Hi; 


e.-y 


^v 


*V 

i ’ 



ROTOR 


2 


(“.“ *'’«')* T Ir, s)' 


(13) 


where Kj^ is the restoring spring constant of the rotor and dependent on the EIj^ 
of the rotor material. 

The complete potential energy for the system now can be expressed as: 


f 

r 


• + [kr " Ir, 

* 


(14) 


whereas the complete rotational kinetic energy for the system is the sum of main 
body, nuiation damper, and rotor components: 




(15) 


and the rather lengthy expression for is presented in Reference 1 (in 

report form). 

The equations of motion for the system can be expressed in terms of the quasi- 
coordinates (<^j, the angle swept out by the pendulous nutation damper 

(0j), and the rotor deflection angles (a^, a^) according to:® 


_d_ 

dt l(9cc>j 


- CO 


dL 
3 dco^ 


+ CO 


dL 


2 dco. 


= L, 


d / 

dt I dco^ 


CO 


dL dL 

1 dco^ ^3 Jco^ ~ ^2 


d / dL 


dL 


dL 


dt [dco^ I “ ^2 doo^ ^ doo. 


= L. 
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(lb) 


Numerical studies with the previously developed nonlinear equations of motion^ 
indicate that whenever the transverse components of main body angular velocities 
are the same order as , the nutation damper motion is characterized by large 
amplitudes and eventually makes contact with the mechanical stops placed on the 
spacecraft to limit displacement to ±20 di^g^ees. Because a stability analysis 
of such a discontinuous system is beyond the scope of the present analysis, it 
vv’as decided to expand Equations (16) under the assumptions that oj. < 1, ^ < 1, 

0J <1, < 1, and ^ < 1. 

The following first order nonlinear equations of motion result: 


+ 0^1 (A-C) + mri (ri + Tq) 01 = (17) 

• 


+ I. 




< + ^0)^ + s]^ sin st “ 

+ s] ^ a^l cos st ► 


J 


y = Ljds) 


- T ^ '1 7 u. 

C^3 + s - mr ^ 01 


+ I, 


cos st + + s] sin st^ = 


L 3 (19) 
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^ *"^2 ‘’ i — n — ■ “ ( 20 ) 

M 


i 

\ 

\ 

Ir 

a + k _ a 

X R z 


s) 2 ] 

a 

z 





1 

i 

) 


^ I 

*■>. 1 

[“1 * ('“'j * 2 s) 

sin 

St + 

“3 * (“2 

+ 2 sj ojjJ cos St 

^ - 0 

( 21 ) 

'i 

i 

i 

Ir 

a + k _ a 

X R X 

* [Kr^Ir, K" 


a 

X 







+ I ^ 

r 

["3 ' (“2 ^ 2 s) 

sin 

St ” 

i-i - (o-j 

+ 2 sj cos st 

^ = 0 

( 22 ) 


where the super barred inertia terms have main body, rotor and natation damper 
components: 


•■4 


A 





B 




For small displacements of the nutation damper, and for m/M <« 1, the system 
center of mass shift due to this motion will be very small and is not included in 
these equations. It is assumed, of course, that the satellite is statically balanced 
in equilibrium when = 0. 


B. Stability Criteria 


Equations (17)— (22) when linearized about the equilibrium motion: cx)^ = fl. 


u), - co^ - d), - a - a - co, - co. - d>. - a - a = 0 

1 3^1 X t 1 3^1 X i'' 

( 

» 

1 would represent a set of nonautonomous differential equations with both constant 

^ and periodic coefficients. The stability of such a system could be analyzed using 

Floquet theor>' similar to the treatment of Mingori.2 The application of Floquet 
^ theory for this problem would necessitate the use of a digital computer to study 

,1 
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the stability point by point by varying different satellite and rotor inertia, spin, 
and damping system parameters in r*. systematic manner. 

In an effort to investigate the stability of this system analytically, the technique 
employed by Flatley*^ in his doctoral dissertation, and applied to the earlier 
problem of Mingori, was applied. 

It is assumed that the effects of the nutation damper and wheoi deflections rela- 
tive to the hub are small perturbations on the nominal rigid body motior^ of the 
torque-free dual spin system. The "zeroth order" equations which represent 
the reference motion, or unperturbed state, can be obtained by assuming that 
terms such as and mr ^ ^ are much smaller than main body inertia terms A, 

B, and C, and also that «A, B, or C. The SAS-A satellite is almost perfectly 
symmetrical about the b^ axis so that A = C; also under the approximations 
mentioned above A = A, B = B. The resulting zeroth order equations corresponding 
to Equations (17), (18), and (19) with Lj = = Lj = 0 are; 



= 0 


(23) 

Aojj 

- 6 ;ja, 3 (B-A) - s 

= 0 

(24) 

Aij 

+ o>j ojj (B - A) + s 

= 0 

(25) 


Equation (23) has the solution = fi, which when substltivced into Equations (24) 
and (25) yields: 


- \ct)^ = 0 


(26) 


a;^ + = 0 


(27) 


where 


A = 



U 


I The solution to Equations (2C) and (27) has the form 


= o)sinX.t 


( 28 ) 


Cl)j = 0)C0S\t 


(29) 


\ 








-•) 



»V'» 


■4 

♦ 


>■ <1 


•I * 



where a; = oj^ (0), the value of at some reference time. 

Thus fo" the unperturbed reference motion of the main body, the angular velocity 
component along the nominal spin axis is a constant while the vector sum of the 
transverse components rotates around b^ with tho nutation frequency K. If s = 0, 
the system reduces to the classical problem of the spinning symmetrical rigid 
body. 

To determine the first order damper and wheel deflection motions all terms 
linear in m and are retained. After substitution of the zeroth order main 
body angular velocities into the equations for the damper and wheel deflections, 
the following equations result: 


+ (k + mr j Tq Q^) = -mr^Paj ^(\ + Q)sin\t] (30) 

"r, * ‘‘r“. * [*‘r * Ir, -n-2s) sin(s -\) t] (31) 

Ir a, + kgi, + Fk^ + + = Ir -n-2s) cos(s -\) t] (32) 

1 L 1 J 1 


Equati )ns (31) and (32) correspond identically to the approximate equations of 
motion for the wheel deflections as derived by Sen and Fleisher,^ 

The steady state solutions to Equations (30)~(32) can be written in the form: 


- Cj sin^t + Cj cos\t 
= Cj sin ( s - A.) t + cos ( s - A.) t 


(33) 

(34) 
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CjSin(s~\)t * cos (s “ X) t 


(36) 


a 


where 

Cj = - mrj rcj(X +fi) mfj (r^ Tj X^)jy/Oj 

Cj ~ mfj rkcjX(X ♦0)/Dj 

C, = Ig^<^\-n-2s)(s->.)k^/Dj 

C, = C, 

^6 " “ ^3 

and 

Dj = [k " ’’ i ^ 


In the absence of the external torques the total angular momentum vector of thc3 
system about an axis passing through the system center of mass (point 0) remains 
time invariant. If the nutation angle, is defined as the angle betveen the b^ 
axis and the total angular momentum vector of the system, H, then 


Kn = Ho cos y 


(36^ 









where = |hJ, Hq = IHq^ h^l. The component Hq may be obtained bv ex- 
panding the following: ^ 



where P,„^ is the linear momentum of the nucaHon damper about an axis through 
0, Vq is given by Equation (2), m is the pendulum end mass, and is the mass 
of the rotor. After exixinsion of Equation (37) and neglecting terms involving 
M '< m 


+ m |(r,2 t * 2r^ r, cos «>,) +i,) 

•> r [-0., sin.#, - O.J (r„ * r, cos#,)j 

• (■'o " * Irj K " 

■ Ir. [^.(-^'jSinst + a., cosst) 


+ 


a 

X 



Cv^S St + 


sin St 


); 


H^cosy (38) 


By differentiating both sides of Equation (38) with respect to time, and substi- 
tuting for Irom the completely nonlinear form of Equation (17), the following 
expression is obtained for H^, : 

“ Isinyy = - mr^ sin^^ 

+ mPr^ ra;^ co^ sin a>^ cos (p^ + (- co^ cos <p^ + sin <p^ )1 
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•f 


+ 





cos St (v^ sin St 



sin st “ 


COS St 



+ 


("j '*) *“.] 

s) -i.] 


+ 




+ 



a a ^ 

% X / 


( 39 ) 


Equation (39) is an exact expression except for the limitation on the magnitude 
of m. If m = = 0, tlie cone angle y is a constant of the motion. 

To analytically express the first order perturbation of the cone angle, the zeroth 
order expressions for main body rates: = a; sin\t, - a; cos \t are 

substituted into Equation (39) to yield: 


H sin 77 = - mr^ 0^^ sinc^j 


mrria;'^[r2 cosM sin0^ - 0 sin\t cos0j j 

(“ sin Xt cos <p^ + cos Xt * sin 
Irj ^ {[cos (s-X)t] [a, (n+s)+a^j 


+ [sin(s-X)t] [o^x 


Substitution for and a^, Equations (34) and (35), would produce sin^ and cos^ 
terms in ^ s - X) t, while substitution of 0^ , Equation (33), would produce sin* 
and cos^ tCims and higher order terms in Xt, after sin 0^ and cos 0^ terms in 
Equation (40) have been expressed in series form. With the use of trigonometiic 
identities, the squared terms may be eliminated, the higher order terms re- 
duced, so that the resulting equation would contain only constants and sinusoidal 
oscillatory terms. 
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Following Flatley,' in the spirit of the method of averaging, the systematic 
change of y is considered to be influenced mainly by the constant terms, with the 
sinusoidal terms contributing only small perturbations of y about some average 
value. 

After redlining only the constant terms in the above averaging process, the 
secula’" motion of the nutation angle can be described by; 


y - ■ 2 H Tiny ["'■'i " Ir, (2s-\ + n> (c, +Cj)j (41) 

Substituting the values of Cj, C^, and from Equations (33)— (35), and approxi- 
mating H sin 7 by the transverse angular momentum of the unperturbed motion. 
Equation (41) becomes: 


y 


H sin 7 
2k^ 


(mij 


21 * (2s +n>2 k^(s -\) 


(42) 


where 


o) - H sin 7/ A 


and 


Q - Th cos 


A 


[^(B- A^ H cos 7 + A I 


sjy^BA 


Eq^iation (42) has the form y - f(7), and equilibrium will be possible whenever 
f( 7 ) =0. Stability at an equilibrium point, 7 ^, requires f(yo) = 0 and 
df/d-) \y <0. Thus, the attitude behavior of this ninth order, nonline i.i system 
is rep resented by the single first order differential equation. Equation (42). 
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Ify = f( 7 ) = F(r) siny, then an equilibrium state will always exist at - 0, 

= n, and, possibly at intermediate nutation angles, if ■ 0. For this ap- 

plica/'ion the nutation angle is bounded by 0 £y £77 . For stability, from considera- 
tion of f'(7o) ^ ^0 ■ necessary that F(0) <0;at>'Q = n, it is nec- 

essary that F(rr )> O; at y^, it is necessary that, F^y^) = OandF'(yQ) < 0. 

C. Special Case- Damper on Spacecraft Only 

For the case where there is only damping present on the main part of the space- 
craft and no damping on the wheel. 


F 


H 

2A^ Dj 


[(mr,r)2(f!+\)2 



(43) 


when k = 0, or = r = 0, or =0, then F = 0 and the system is unstable for 
all values of >q. The conditions that, for stability at y^ = 0, k > 0, ^ > 0, > 0 

are part of the necessary and sufficient stability criteria previously derived by 
the methods of Routh-Hurwitz for the case of spacecraft damping only.^ Assuming 
these conditions are satisfied, the sign of F is clearly determined by the sign of 
Therefore, for stability at y^ = 0, - \ > 0, or ^ > 0, or 




CB-A)H^^ + AI,^s 
AB 


> 0 


(44) 


Condition (44) can be shovm to be equivalent to: 


B + - A > 0 (45) 

within the assumptions previously made regarding A, A, and B, B. Condition (45) 
is equivalent to inequality (16) of Reference 1, resulting from the Routh-Hurwitz 
analysis. 
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It is apparent that from consideration of Equation (43), the necessary and suf 
ficient stability criterion, (18), of Reference 1, repeated here: 


0 m _ mP 

- ^ + K > 

1 


M 


B-AM 




( 46 ) 


is ^ apparent. This can possibly be exi^lained by the fact that m/M and (for 
small or moderate 1) are higher order terms when compared with. A, B, or 
Ij^ sy/o, as in the present analysis. In practical application, condition (46) limits 
the maximum height of the damper plane for a given set of main body inertias, 
nutation damping system parameters, and ratio of s/Q. 


From consideration of Ineq. (44), the line = 0 or = (A- B)/Ir 2 » divides 
the initial state plane into the regions \ 0, <0. The region: > 0 is sta- 

bilizing for values of 7 ^ = 0 , whereas the region; < 0 describes stability for 
- 77. Furthermore, the line: BO^^I^^s = Hcos 7 = 0 divides the (Q^, s) 
plane into regions describing y^ = 0, and 7 q = Superimposing the two results, 
the following stabiliU- diagram can be constructed in the (0^, s) plane (Figure 3), 
where the shaded regions indicate stability. 

As Flatley indicates^ a useful alternate form of this stability diagram is obtained 
by letting: q = B A ~ B A, andq' = Ir^A ~ The boundary line, 

s/Qq = (A- B)Ir2 becomes 1 = q + q' s/Q^ and the boundary line; BOq + Ir^ s = 0 
becomes q+ q' = 0. Th'=i transformed stability diagram appears in Fig- 

ure 4, where the physically meaningful values of q are restricted to q' < q < 2. 
Unless the momentum wheel is deliberately designed to spin opposite to the 
sense of the main body angular momentum vector, the physically meaningful 
values of are to the right of the q axis. Figure 4 is identical to a result 
obtained by Mingori^ using some of the Routhian stability boundaries. 

A critical value of rotor spin*^ occurs when only one stable equilibrium nutation 
angle exists and = 77 is unstable. This value may be obtained by setting ^ = 0 
when 7 q = 77, yielding 


"cH, = (B-A)H„/aI„^ (47) 

When s > s^j.. £.t y ^ =77, then \ becomes positive, and thus F(t 7) < O; therefore, 

7q = 77 becomes a point of unstable equilibrium, i.e. a dual-spin satellite which 
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STABLE 


UNSTABLE 


UNSTABLE 


BOo+Ir s-0 


Figure 3. Stability Diagram With Damping Only on the Spacecraft 


has inverted itself from its nominal pointing direction, will, in the presence of 
main body damping, tend to restore itself to the "upright” position. 

D. Special Case-Damper on Momentum Wheel Only 


If there is only damping present on the momentum wheel, 


F = 


(2s-\ + fi)^k„(s-X)] 


when = 0, then F = 0, and the system is unstable for all values of . For 
this special case, if ^ 0 and there is damping on the wheel, the sign of F is 
directly dependent on the sign of { s-X ). An initial nutation angle of zero is 
stable if (s - Xq) is negative, while an initial nutation angle of 180*’ is stable if 
(s - \q) is positive. Therefore, the line " s = 0, or s/n^ = (B - AV^(A - 






Figure 4. Stability Diagram in q, q* Space With Damping Only on the Spacecraft 


divides the s, plane into the two possible regions: “ s >0 and - s < 0. 

Recalling that the line s/0^ = - b/Ir^ sefjarates the 0^, s plane into the two re- 
gions correspondirg to - 0 and >q = tt respectively, a stability diagram may 
be constructed (Figiare 5). 

Again tliis diagram may be reconsiructed in terms of the nondimens ional param- 
eters <1 and q', as before. The physically realistic values of q are given by: 
q' ' q < 2. Figure 6 illustrates the alternate form of Figure 5 for the special 
case of wheel damping only. 

E. Discussion of General Case 


The general expression for time rate of change of nutation angle, (42), may be 
expressed as: 






y 


2H sin 7 


D. 


(2s-X + n)k^(s-\) 

+ IT 


(49) 
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UNSTABLE 



Fig-jre 6. Stability Diagram in q, q' 


Space With Damping Only on the Wheel 


after su’.jtituting 1/A^ - sin^y. From this, necessary conditions for 

stability at the equilibrium positions ofy^ = 0, and 7 q = ^ may be obtained by 
examining the sign of F(7o)» intermediate position, 7 q, by examining 

the sign of F'(7 q) as explained earlier. 

Likins® has previously derived attitude stability criteria for a dual spin space- 
craft in terms of the average energy dissipation rates on both the wheel and the 
despun portion using an approximate energy sink analysis. In an effort to com- 
pare the form of Equation (49) with Likins' criteria the average energy dissipa- 
tion rates were evaluated for the SAS-A dual spin system. 

The average energy dissipation rate due to the action of the pendulous nutation 
damper is: 


a ve 


= “ k0j' 


■ ■ 


(50) 


f’rom Equation (33), 


01^ = cosUt - Cj sin 2U +C 2 ^ sin2 \t) (51) 


In order to obtain , Equation (51) is averaged using the same procedure 

* s s a vt _ 

as before, by wTiting terms involving cos^ and sin^ in terms of constant terms 
and sinusoidal terms, and retaining only the constant terms in the averaging 
process. There results, 





(52) 


After substitution of (52) into (50), 


T = r)2o>2 (\ + Q)k/2Dj 

ve ' 


(53) 
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Similarly for the rotor, the average energy dissipation rate is the sum of the two 
components of average dissipation rates. 


= T. 


■ V* 


+ T. 


a ve 


where 


(54) 



■ V • 




■ ve 



■ ve 




■ ve 


Following the same procedure as above, it can be shown that. 


i 3 


+ 


■ ■ 


■ ■ 


= (s-\)» (C3’ 


(55) 


■ V e 


and, therefore. 


= - (s - (\ -0- 2 s)2 


■ ve 


(56) 


By comparison of Equations (49), (53), and (56), it can be seen that yean be 

e • 

written in terms of T„ and , as follows: 

"eve *^eve 


y 


1 

H siny 



(57) 


or 


d 

dT(Hcosr).,. = 


■ ve 


\ “ S 


(58) 
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The terms within the bracket of Equation (57) are identical in form to those ap- 
pearing in Likins' energy sink analysis. ^ As Likins indicates the terms Td,v«A 
and -s)are the "first approximation" dissipative torques about the 

spin axis of the main body and spin axis of the rotor respectively, which is ap- 
parent from consideration of Equation (58). The implication here is that if the 
averaging technique of Flatley'^ were carried to a higher order, (e.g. first order 
approximation for main satellite angular rates, and second order approximations 
for dami^er and wheel deflections), additional terms may result in F(y), yielding 
additional, more complete stability Information. 

in. NUMERICAL RESULTS 
A . Calculation of Stability Oiteria for SAS-A 

For stability, from consideration of f' < 0, at^Q = 0* it is necessary that 
F(0) <0. F(0) was numerically evaluated using representative SAS-A param- 
eters as listed below: 


27.0 kg-m ^ 

k 

= 7.0 X 10"* nt-m-sec/rad. 

28.54 kg-m^ 

K 

= 6.10 X 10'* nt-m/rad. 

0.44 m. 

n 

= 8.72 X 10'^ rad/sec. 

0.019 m. 

m 

= 0.2158 kg. 

0.1833 rn. 


= 71.6107 nt-m/rad. 

0.011519 kg-m^ 


= 0.006778 nt-m-sec/rad. 

0.0057595 kg-m2 




*'1 

Four different values of s were considered: 0, 0.1 nominal, 0.5 nominal and the 
nominal value of 209.4 rad/sec. 

We write F(y ) as 

F(r) = +Q,] (58) 


•Determined experimentally at The Applied Physics Laboratory, Johns Hopkins University. 
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where 


■ ■ [(">f 1 1')^ •‘'^jy 


= [2I,» (2s-X + n>»k^( 


The results are summarized in the following table. 



s/Oq = 


0 

2401.3 

12,006.88 

24,013.7615 


Qi (0) 

Qj <0) 

[2A^ h] • F(0) 

-2.448 X 10“’ 

-2.949 X 10"*® 

-2.448 X 10*’ 

-1.8379 X 10"* 

3.2185 X 10"* 

-1.516 X 10"® 

-1.3167 X 10"* 

4.015 X lO""* 

-0.9152 X lO"-* 

-0.4092 

3.195 X 10"^ 

-0.406 


The system is stable for all four values of q , but the margin of stability tor tne 
two intermediate values of q' is not as groat as at the nominal ofxjrating con- 
dition where the stability is insured by a factor of 128. This calculation does not 
guarantee the system stability during the actual spin-up process, but does give 
an indication of stability if, for some reason, the wheel does not reach its nominal 
operating angular velocity, s. (To study the problem of stability during spin-up, 
an equation for the motor torque would have to be included in the formulation 
with s now treated as a function of time instead of as a constant.) The sign of 
Oj (0), in the numerical example, verifies that only at very low values of q', is 
the effect of wheel damping. stabilizing. The "X" appearing on the stability charts. 
Figures 4 and 6, refers to the SAS-A satellite under nominal operating conditions. 

B. Calculation of 

Using the SAS-A parameters listed in III. A., the critical value of wheel soin was 
calculated from Equation (47) to be ^crit~ 26.35 rad/sec. Since the nominal 
value for s under normal SAS-A operating conditions is 209.4 rad/sec, s > . . ., 















iind the nutation angle of 180’ (satellite inverted) becomes a point of unstable 
equilibrium, if the effects of wheel danjping nre neglecbKl. This moans that should 
SAS-A reach m. ir vei tod posltica, gi /eii i sit. all iM}rt».i b ilion, it would tend to 
restore Itself to the ’ upright” iwsition. (In this argument the effect of such dis- 
continuities as nutation damper stops have not been included, but could be de- 
termined using a digital computer simulation of the nonlinear equations of motion 
with a suitable model for the stops.) 

C. Time Consta nt A ssociated With Nutation Angle Decay for 

The time constant associated with nutation angle decay is given by the reciprocal 
of the slope of f(y) a^ ;ne equilibrium point. Equation (42) has the form r = f(7) 
and at an initial point of equilibrium, f (yg) ~ If we let y = y^ + x where 

X is small when compared with y^, then, y “ x. and f(y) ^ (^ o) 

Thus, X % xf'(>j,) and dx x = f'(yQ)dt. This differential equation has the 
exponential solution of the form: x = x(0) = x(0) e* where t is the 

exixmential time constant, and t = - 1' f ' (yg)* 

For SAS-A we are interested in the time constant associated with nutation angle 
decay for an initial equilibrium nutation angle of zero degrees. At yg = 0, 
f' (0) = F(0), and r = - 1 F(0) . Using SAS-A parameters it was found that 
r = 22.3 minutes, which compares favorably with previous numerical integration 
results not including wheel damping.^ 


rV\ CONCLUDING REMARKS 

Asa result of the present stability analysis and numerical results the following 

conclusions can be made: 

1. The stability criteria obUiined for this system by means of ave^'Jjging'^ the 
"first order" perturbations on the system nutation angle can be reduced to 
the criteria emanatirz from Likins’ energy sink analysis. The implication 
is that higher order - -.ns would have to be also included in the averaging 
process to yield mo • r /.mplete stability Information. 

2. In the absence of wheel damping, some (but not all) of the previously derived^ 
Routh-Hurwitz necessary and sufficient conditions can be obtained from the 
present analysis. 

3. In the absence of wheel damping, the nominal SAS-A rotor spin rate is such 
th: t a nutation angle of 180 degrees (inversion) is a position of unstable 
equilibrium. 
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4. A numerical evaluation of the analytical stability criterion developed here, 
using SAS-A satellite parameters and measured wheel damping data, indi- 
cates that cystem stability is insured by a factor of 128 under the normal 
operating conditions. The system is also stable at smaller values of rotor 
spin rate, although the margin of stability is somewhat reduced. 

5. The time constant associated with nutation angle decay under me normal 
SAS-A operating conditions is 22. .1 minutes, and ^vs not appreciably degraded 
by the measured energy dissipation in the wheel. 

The nonlinear equations of motion derived here will be solved by numerical in- 
tegration techniques in order to 1) verify the results of the present analytical 
study, 2) consider further stability problems during momentum wheel spin-up, 
and 3) evaluate the effects of external perturbations (such as aerodynamic and 
gravity gradient torques) on the results reported here. The results of these 
numerical studies will be reported subsequently. 
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APPENDED A 
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The following two figures show the results of wheel damping tests performed at 
the Applied Physics Laboratory by Mr. R. E, Fischell and Mr. W. Radford. The 
values of and were determined for both air and vacuum conditions. With 
the wheel initially stationary (Figure A.l) and then initially spinning (Figure A. 2) 
and with the wheel symmetry axis vertical, an initial impulse v/as applied to the 
outer rim of the wheel. Instrumentation recorded the vertical displacement of 
the wheel as a function of time. It can be seen that the transient responses shown 
in the Figures are, for all practical purposes, characteristic of a second order 
linear system. Values of K^and selected for numerical examples in this 
paper represent the most pessimistic case (worst destabilizing effect). 

The wheel deflection Equations (21) and (22) are of the form; 

; 4 ^; ■ 

^ * [Kr+Ir, = f(t) (Al) 

which have a ti'ansient solution, 


a - e [^sin t + (//q)J (A2) 

where O q and are determined from the boundary conditions. The time constant, 

T, is the value of t when kj^ ^/2Ip = 1. Thus 

kR = 21^/r (A3) 

The damped natural frequency, 0^, is related to the system parameters according 

K 

to: 







Kr"Ir, K" 


VIr J " 
\ V 



— J 



1/2 


(A4) 


1/2 

and \ is the undamped natural frequency. 







